Phase space analysis of quintessence fields trapped in a Randall-Sundrum Braneworld: 
anisotropic Bianchi I brane with a Positive Dark Radiation term 
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In this paper we investigate, from the dynamical systems perspective, the evolution of an scalar 
field with arbitrary potential trapped in a Randall-Sundrum's Braneworld of type 2. We consider 
an homogeneous but anisotropic Bianchi I (BI) brane filled also with a perfect fluid. We also 
consider the effect of the projection of the five-dimensional Weyl tensor onto the three-brane in 
the form of a positive Dark Radiation term. Using the center manifold theory we obtain sufficient 
conditions for the asymptotic stability of de Sitter solution with standard 4D behavior. We also 
prove that there are not late time de Sitter attractors with 5D-modifications since they are always 
saddle-like. This fact correlates with a transient primordial inflation. We present here sufficient 
conditions on the potential for the stability of the scalar field-matter scaling solution, the scalar 
field-dominated solution, and the scalar field-dark radiation scaling solution. We illustrate our 
analytical findings using a simple /-deviser as a toy model. All these results are generalizations of 
our previous results obtained for FRW branes. 
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I. INTRODUCTION 

Randall-Sundrum braneworlds were first proposed in 
[l|, HJ. In these references was proved that for non- 
factorizable geometries in five dimensions there exists a 
single massless bound state confined in a domain wall 
or three-brane. This bound state is the zero mode of 
the Kaluza-Klein dimensional reduction and corresponds 
to the four-dimensional graviton. The Randall-Sundrum 
brane of type 2 model, as an alternative mechanism to the 
Kaluza-Klein compactifications [l[ , have been intensively 
studied in the last years, among other reasons, because 
its appreciable cosmological impact in the inflationary 
scenario 0-(5|. The setup of the model start with the 
particles of the standard model confined in a four dimen- 
sional hypersurface with positive tension embedded in a 
5-dimcnsional bulk with negative cosmological constant. 
It is well-known that the cosmological field equations on 
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the brane are essentially different from the standard 4- 
dimensional cosmology pHs)]- 

Friedmann-Robertson branes with an scalar field 
trapped on it have been investigated widely in the lit- 
erature 

§£1. 

In [14] was shown that the potential 
V oc cosech 2 (A<j>) leads to scaling solution in a RS2 sce- 
nario. The dynamics of a scalar field with constant and 
exponential potentials was investigated by [TBI ]. These 
results were extended to a wider class of self-interaction 
potential in [l6j using a method proposed in [l7j sup- 
porting the idea that this scenario modifies gravity only 
at very high energy/short scales (UV modifications only) 
having an appreciable impact on primordial inflation but 
does not affecting the late-time dynamics of the Universe 
unless if the energy density of the matter trapped in the 
brane increase at late times [l8[ . In the reference [l9| we 
investigate a scalar field with arbitrary potential trapped 
in a RS2 Braneworld. There we present sufficient condi- 
tions for the asymptotic stability of de Sitter solution 
and for the stability of scaling solutions as well as for the 
stability of the scalar-field dominated solution extending 
the results in [2(| to the higher dimensional framework. 
We prove the non-existence of late time attractors with 
5D-modifications -a fact that correlates with a transient 
primordial inflation. Finally, as an example we studied a 
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scalar field with exponential potential, 
V(4>) = V e~ x * + A 



(1) 




20( . and in the 
In both cases 



confined in the brane. We proved that for x < the 
de Sitter solution is asymptotically stable. However, for 
X > we have proved the de Sitter solution is unstable 
(of saddle type). The exponential potential ([T]) have been 
widely investigated in the literature. It was studied for 
quintessence models in [21] were it is considered a nega- 
tive cosmological constant A. In our case of interest we 
assume A > 0, to avoid dealing with negative values of 
y. However we can apply our procedure by permitting 
negative values for y for the case A < 0. The dark en- 
ergy models with exponential potential and negative cos- 
mological constant were baptized as Quinstant Cosmolo- 
gies. They were investigated in (22j using an alternative 
compactification scheme. The asymptotic properties of a 
cosmological model with a scalar field with exponential 
potential have been investigated in the context of the 
General Relativity by the authors of 
context of the RS2 braneworlds by [l5 
it was studied of the pure exponential potential (A = 0). 
Potentials of exponential orders at infinity were studied 
in the context of Scalar-tensor theories and conformal 
F(R) theories by the authors of [H-Hl]. 

Here we give a step further by considering the natural 
generalization of FRW, i.e., Bianchi I cosmologies. The 
results in fl9j concerning the exponential potential apply 
here too. 

Bianchi I models are the minimal extension of FRW 
metric to the anisotropic framework. Homogeneous but 
anisotropic geometries are well-known [27l r28j . Bianchi 
I, Bianchi III, and Kantowski-Sachs can be a very good 
representation for the homogeneous but anisotropic uni- 
verse. They were investigated in the framework of f{R) 
cosmology from both numerical and analytical view- 
point also incorporating the matter content (see [29| and 
the references therein). The evolution of cosmological 
braneworld models were investigated, for instance, in 
[30l - l35j . In [32| it is presented a systematic analysis of 
FRW, Bianchi I and Bianchi V metrics in these scenar- 
ios. There it is discussed the changes in the structure of 
the phase space with respect the general-relativistic case. 
In [36| it is studied the dynamics of a BI brane in the 
presence of inflationary scalar field and it is showed that 
the high energy effects from extra dimensional gravity 
removes the anisotropic behavior near the initial singu- 
larity which is found in general relativity. However, if an 
anisotropic stresses are included the model behavior in 
the vicinity of the initial singularity changes. In [37] a 
BI brane-world model with a pure magnetic field and a 
perfect fluid with a linear barotropic 7-law equation of 
state was studied being found that if 7 < | 1 the past 
attractor corresponds to a critical point with non trivial 



1 If 7 > I the previous fixed point is not the only past attractor 



magnetic field, which is also a local source, bringing the 
anisotropic behavior back in the initial singularity. 

In this paper we are interested in investigate the evolu- 
tion of an scalar field with arbitrary potential trapped in 
a Randall-Sundrum's Braneworld of type 2. We consider 
an homogeneous but anisotropic Bianchi I (BI) brane 
filled with a mixture of an scalar field with arbitrary po- 
tential and a perfect fluid with equation of state (EoS) 
parameter u = 7 — 1 with 1 < 7 < 2. To the potential 
treatment we use the method developed in [l7|. We ex- 
tend previous results in [IH, [l|| [l9|, [23|, [32j by considering 
scalar fields with arbitrary potentials in an anisotropic 
background and considering the effect of the projection 
of the five-dimensional Weyl tensor onto the three-brane 
in the form of a positive Dark Radiation term. In or- 
der to illustrate our analytical results we consider a toy 
model using the simple /-deviser given by f(s) = s + a. 
The cosmological implications of such a model are also 
discussed. 



II. BIANCHI COSMOLOGY IN THE BRANE 

The set up is as follows. Using the Gauss-Codacci 
equations, relating the four and five-dimensional space- 
times, we obtain the modified Einstein equations on the 
brane [SI, El 2 : 



G a b 



-A 4 g Q 6 + k T a 



"(5) 



S a b &ab 



(2) 



where g ao is the four-dimensional metric on the brane 
and G a b is the Einstein tensor, k is the four-dimensional 
gravitational constant, and A4 is the cosmological con- 
stant induced in the brane. S a b are quadratic correc- 
tions in the matter variables. Finally, £ a b are corrections 
coming from the extra dimension. More precisely, £ ao 
are the components of the electric part of the Weyl ten- 
sor of the bulk (see [3^| and the review [Ilj]). Following 
[231 . 13TI H2| , from the energy-momentum tensor conser- 
vation equations (V a T h a = 0) and equations © we get a 
constraint on S a b arid E a b'- 



7°" (y£ab ~ K^S a b 







(3) 



In general we can decompose £ ao with respect to a chosen 
4- velocity field w a [42| as: 



£ab 



U(u a u b + -h ab ) + V ab + 2u (a Q h ) 



(4) 



possibility. It was argued, in both cases, that chaotic behavior is 
not possible [37l[3i|. 
2 In this section, we follow \'2'-i |3ll in order to give a brief outline 
of the considerations that lead to the effective Einstein equations 
for Bianchi I models 1191 14H . 
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where 

V {ah) = Tab, V% = 0, V ab U b = 0, Q a U a = (5) 

here the scalar component IA is referred as the dark en- 
ergy density due to it has the same form as the energy- 
momentum tensor of a radiation perfect fluid. Q a is an 
spatial vector that corresponds to an effective nonlocal 
energy flux on the brane and V a b is an spatial, symme- 
try and trace-free tensor which is an effective non local 
anisotropic stress. Another important point is that the 9 
independent component in the trace- free £ n h are reduced 
to 5 degree of freedom by equation © [4l],|4j|], i.e., the 
constraint equation Q provides evolution equations for 
hi and Q a , but not for V a b- 

Taking into account the effective Einsteins equations 
©, the consequence of having a Bianchi I model 3 on the 
brane is [36j : 

Qa = (6) 

but we do not get any restriction on V a b [H, Hl[ . Since 
there is no way of fixing the dynamics of this tensor we 
will study the particular case in which: 

Vab = 0, (7) 

this condition, toghether with © and ([3]), implies: 

D a U = 0^U =U(t) (8) 

Using the above conditions over Q a and V a b (00 , set- 
ting the effective cosmological constant in the brane to 
zero, i.e., A4 = 4 , the effective Einstein equations @ 
for Bianchi I models (which have zero 3-curvature, i.e., 
_R( 3 ) = 0) become: 

" 2 = rM'+§)+r 2 +? < 9 > 



*=-K 1+ f)(* 2+ ^)-f-' 2 < 10 » 

a = -3Ha (If) 

U = -AMU (12) 

p m + 3H{ Pm + Pm ) = (13) 

4> - 3H<j) + d^V = 0, (14) 

3 the metric in the brane is given by ds 2 = —a 2 (t)dt 2 + 



4 The induced cosmological constant in the brane can be set to 
A4 = by fine tuning the negative cosmological constant of the 
AdS§ with the positive brane tension A > |4ll |44H . 



where a is the shear term. 

It is also convenient to relate p and p by 

P= (7 -!)/>■ (15) 

The Dark Radiation term in ([ MTU)) evolves as U (t) = 
a ^ 4 [42| . where C is a constant parameter. An impor- 
tant issue comes from the sign that can take C . From 
the point of view of the brane (brane base formalism) C 
is just an integrations constant and can take any sign. 
However it was shown, from the bulk base formalism, 
that for all possible homogeneous and isotropic solutions 
on the brane, the bulk spacetime is Schwarzschild-AdS 
being possible to identify C with the mass of a bulk 
black hole (p = and use this to constrain to be 

positive [43|, |4j| |46|. For anisotropic models such con- 
straints do not exist, i.e., C, and therefore U, can take 
any sign [H, |3l|, HtJ. In this first approach for in- 
vestigating arbitrary potentials, we remain in the case 
C > =^ U > 0. The case of hi < is also wor- 
thy of investigation and deserves a careful analysis. For 
investigating these classes of models, i.e., with U < 0, 
it is requiered to consider a new set of variables than 
in the present paper (more precisely, normalizing with 

D = J H 2 — ZLL instead of H and reemplacing the ana- 
logue to the variable y but the new one Q — |y) 5 . Since 
the resulting system have its own subtetlies and also have 
a rich cosmological behavior, these results would enlarge 
the present report, so we prefer to address the point in a 
forthcomming paper [48|. 

III. DYNAMICAL SYSTEMS ANALYSIS FOR 
BIANCHI I BRANE WITH POSITIVE DARK 
RADIATION TERM 

A. Dynamical variables and dynamical system 

In this section we investigate the cosmological model 
(|9|)- ([T4| from the dynamical systems perspective. First 
we need to recast the cosmological equations (li?|)-([T4"f as 
an autonomous system of first order ordinary differential 
equations [25|, [5^ |49|, [5(| ■ For this purpose we introduce 
the normalized variables 

n m = — 2 = n u = JM ,, (i6) 

the additional dynamical (non-compact) variable, s, 
given by 

s=-d <f> \nV(<f>). (17) 



J This allows for the construction of a phase space compact in all 
the variables but s defined in eq. I|17|l. that runs, in principle, 
over all the real values. 
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which is a function of the scalar field and the new time 
variable t = lna(f). 

For the scalar pot ential treatment, we proceed follow- 
ing the reference [l7l [l9| . Let be defined the scalar func- 
tion 



1, 



V V 



(18) 



Since T is a function of the scalar field T(4>) (see definition 
tH])), also is the variable s = S((j>). Assuming that the 
inverse of S exists, we have </> = S (s). Thus, one can 
obtain de relation T — T(S~ 1 (s)) and finally the scalar 
field potential can be parameterized by a function /(s). 

Using the Friedmann equation (fT0|) we obtain the fol- 
lowing relation between the variables (fT6|) 

x 2 + y + a m + n x + n 2 a + n v = i (19) 

The restriction (|19|) allows to forget about one of the 
dynamical variables, e.g., to obtaining a reduced dy- 
namical system. From the condition < flu < 1 we have 
the following inequality 

< x 2 + y + fl m + fix + ttl < 1 (20) 

Using the variables ([T6|) . the field equations (|9|)-(fT4|) 
and the new time variable we obtain the following 
autonomous system of ordinary differential equations 
(ODE): 



sy + x 3 + [fll - 2y 



3l 



2 n, 



3((7-2)n m -2y)n> 
x 2 + y + fl m 



• 4fi A - 1 



y' = y (2x 2 - V6sx + 2f2 2 - Ay + 3jn m - Ail, 
2 (4a; 2 - 2y + (3 7 - 2)fi ro ) tt x 



x 2 + y + f2 m I 
n' m = il m (2x 2 - Ay - 3 7 + (37 - 4)fi m + 4) + 

n m (8x 2 -A y + (6 7 - A)n m ) n x 



8x 2 



x 2 + y + fi m 
Ay + (67 - A)fl m ) n 



2 £7 777, 



(21) 



(22) 



(23) 



y + ^r, 



9(2 o . 3((7-2)x 2 + j/ 7 ) , „ , „ 
2 a; — 2y — 37 H — + 2 | S2 A + 



y + 



(2f2 2 + (3 7 - 4)fi m ) f2> 



ft = n! 



.x 2 - 2y 



(4x 2 - 2y + (3 7 - 2)Q m ) f2 A » CT 
x 2 + y + VL m 



(24) 



(25) 



s' = -\/6s 2 xf(s) (26) 

Using the Friedmann equation (fT9"]l we obtain the fol- 
lowing useful relationship 

Pt _ fix 
A f2 m + a; 2 + y 



(27) 



From (|27|) follows that the region Q m + a: 2 + y = cor- 
responds to cosmological solutions where pr ^S> A (corre- 
sponding to the formal limit A —> 0). Therefore, they are 
associate to high energy regions, i.e., to cosmological so- 
lutions in a neighborhood of the initial singularity 6 . Due 
to its classic nature, our model is not appropriate to de- 
scribing the dynamics near the initial singularity, where 
quantum effects appear. However, from the mathemati- 
cal viewpoint, these region (fi m + x 2 + y = 0) is reached 
asymptotically. In fact, as some numerical integrations 
corroborate, there exists an open set of orbits in the phase 
interior that tends to the boundary fi m + x 2 + y = as 
t — > —00. Therefore, for mathematical motivations it 
is common to attach the boundary Q m + a; 2 + y = to 
the phase space 7 . On the other hand the points with 
(fi A = 0) are associated to the standard 4D behavior 
(Pt < A or A -> 00) and corresponds to the low energy 
regime. 



From definition (|16l) and from the restriction (|19p , and 
taking into account the previous statements, it is enough 
to investigate to the flow of (|2"Tj) - ([2l))) defined in the phase 
space 



9 = {(x,y,Q m ,ilxM :0< 



ft 



+n x + < 1, -1 < a; < 1,0 < y < i, 

< n m < 1,0 < A < 1, -1 < f2 CT < 1} x {s G K} .(28) 
B. Phase space analysis 



The critical points of the system (j2T|) - (|26|) are summa- 
rized in Table U where s* corresponding f(s*) = 0. All 
the points showed in table Q] but P^, P^ 1 for x = 0, and 
P10 when y = always satisfy the condition £l m +x 2 +y 7^ 
0. The points where this condition do not holds will be 
excluded from our analysis, since they represents solu- 
tions near the initial singularity. As we commented in 
the last section, our model is not applicable near this 
singularity. 

It will be helpful to have the important observational 
parameters in terms of the state variables. The dimen- 
sionless scalar field energy density parameter Q,^ = gjjr, 
the equation of state parameter w<j, = — , and decelera- 

jr) . They read: 



tion parameter q = — 1 



f^0 = x 2 + y, 



^4, = 



y 



-l + 3x 2 



37 



3 (2x 2 + 7 n m ) 

x 2 +y + Q m 



(29) 



(30) 



6 See the references |24l , l5l| for a classical treatment of cosmological 
solutions near the initial singularity in FRW cosmologies. 

7 We submitt the reader to the end of section [lH Bl for a discussion 
of how to dealing with the singular points at the hypersurface 
n m + x 2 + y = 0. 



5 



TABLE I: Existence conditions for the critical point of the system (|21|) - (|26[ ). We use the notations s c for an arbitrary real 
s-value and s* for an s-value such that f(s*) = 0. 



Pi 


X 


y 








s 


Existence 



Pi 








i 












s c G R 


1 2 


±1 



















always 


3 


±1 
















s* 


always 


p. 

" 4 


s* 


3(7-2)7 


1 — 3 7, 


n 
u 


n 
u 




s 


*2 \ q 


Pr, 


s 

vli 


.* 2 
1 — 













s* 


s* 2 < 6 


P ± 


vs 













6 


s* 


s* 2 > 6 


Pi 


VI 

s* 


4 

3s* 2 













s* 


s* 2 > 2 


P ± 


x-G [-1,1] 











±Vi- 


a; 2 







P ± 


xe [-1,1] 











±Vi- 


X 2 


s* 




Pro 





i-n A 





fix 6 [0,1[ 










G [0, 1 


Pll 











1 







s c g R 


always 



In the table U are displayed the existence conditions for 
the critical point of the system (|21 [) -(|26 [1 whereas in the 
table |n] are shown the values of the cosmological param- 
eters uty, f2^ and q. 

Now, let us comment on the stability and physical in- 
terpretation of the critical points/curves. 

The point P± represents a matter-dominated solution 
(Q m = 1). Although it is non- hyperbolic, it behaves like 
a saddle point in the phase space of the RS model, since 
they have both nonempty stable and unstable manifolds 
(see the table HH]). 

The critical points P^ are solutions dominated by the 
kinetic energy of the scalar field and they represent so- 
lutions with an standard behavior = 0). This crit- 
ical points are nonhyperbolic. However, they behave as 
saddle-like points in the phase space because of the insta- 
bility in the eigendirection associated with two positive 
eigenvalues and the stability of an eigendirection associ- 
ated to a negative eigenvalue. 

The critical point are solutions dominated by 
the kinetic energy of the scalar field. They are non- 
hyperbolic, however they behave as saddle points since 
they have both nonempty stable and unstable manifolds 
(see the table ITTT) ) . Thus, they represent transient states 
in the evolution of the universe. 

The critical point P4 is nonhyperbolic for 7 = 2,4/3, 
s* = ±\/37 and f'(s*) = corresponding to scalar field- 
matter scaling solutions (f2^ ~ £l m )- 

P4 is a stable node in the cases < 7 < |,s* < 

-V37./V) < or § < 7 < i-jfy < < 



-V37,/'(0 < 0, or < 7 < Is* > V37, /'(«*) > 0, 
or § < 7 < iV^< s* < *fi^J>(s*) > 0. 

It is a spiral stable point for either |<7<|,s*< 

-7§fe,/'(**) < or § < 7 < l,s* > > 
0. 

In summary, it is stable for either < 7 < |,s* < 

-V37, /'(**) < or < 7 < f,s* > V37, /'(«*) > 0. 
Otherwise, it is a saddle point. 

The critical point P5 represents a scalar-field dom- 
inated solution (51^ = 1) that is not hyperbolic for 
s* e {0,±V6,±x/3t",2} or f'(s*) = 0. 

In the hyperbolic case, P5 is a stable node for either 
< 7 < 1,-^37 < s* < 0, /'(«*) < or § < 7 < 
2,-2 < s* < 0,/'(s*) < or < 7 < |,0 < s* < 
y/5y,f'(s*) >0or|< 7 <2,0<s*< 2J'(s*) > 0; 
otherwise, it is a saddle point. 

The critical points P^ 1 are nonhyperbolic and corre- 
spond to scalar field-anisotropic scaling solutions (0^ ~ 
Via). For <C 1 the anisotropic term in the Fried- 
mann equation (ITU1) dominates the cosmological dynamic 
(fi CT = ±1). However, they behave as saddle-like points 
since they have both nonempty stable and unstable man- 
ifolds (see the table IIII[) . 

The critical point P7 is nonhyperbolic for 7 = 4/3, 
s* = ±2 and f'(s*) — and corresponds to scalar field- 
dark radiation scaling solutions (f2^ ~ &u)- 

It is a stable node for either | < 7 < 2, — -7= < s* < 
-2,f'(s*) < or § < 7 < 2,2 < s* < /'(«*) > 0. 
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TABLE II: Some basic observables S7<£ and q for the critical point of the system (|21|I -(I26 | ). We use the notations s c for an 
arbitrary real s- value and s* for an s- value such that /(s*) = 0. 



Pi ^<jy 




q 




P, 
1 1 


1 1 Tl H pfi Tl pH 

U.lH_lt_.±l±±t_.LJ. 


o 


37 1 
2 1 


1 2 


1 


i 


2 


P± 


i 
i 


i 
i 


o 
I 


p, 


7 1 


3t 

' - 


37 1 

2 ~~ 1 


J 5 


3 1 




2 i 




1 


(> 




Pr 


1 
3 


4 
>- ' - 


4i-x 


P ± 


1 


a; 2 


3a; 2 - 1 


P± 
-"9 


1 


x 2 


3a: 2 - 1 


Pio 


-1 




-1 




undefined 





undefined 



It is a spiral point for either | < 7 < 2, s* < 
-$g,f (**) < or § < 7 < 2, S * > -8.,/' (a*) > 0. 

In summary, it is an attractor for either | < 7 < 
2, s* < -2, /' (s*) < or § < 7 < 2, s* > 2, /' (s*) > 0; 
it is a saddle otherwise. 

The circles of critical points and P^ are nonhyper- 
bolic and corresponding to scalar field-anisotropic scaling 
solutions (f2^ ~ r2 CT ). Both solutions represents transient 
states in the evolution of the universe. When x — > the 
anisotropic term in the Friedmann equation (jlOl) dom- 
inates the cosmological dynamics (fi CT = ±1). In this 
limit the corresponding cosmological are in a vicinity of 
the initial singularity. 

The line of critical points Pio is nonhyperbolic. They 
represent solutions with 5D-corrections for Q\ =/= 
whereas for £1> = it represents an standard 4D cos- 
mological solution. From the relationship between y and 
Q\ follows that this solution is dominated by the poten- 
tial energy of the scalar field pr = V(4>); that is, it is 
de Sitter-like solution (w^ = — 1). In this case the Fried- 
mann equation can be expressed as 

3H2 = V { 1 + Jx) ^ 

In the early universe, where A<7, the expansion rate 
of the universe for the RS model differs from the general 
relativity predictions 



Hrs 
Hgr 




Due the importance of de Sitter solutions in the cosmo- 
logical context, in the section HVl we calculate explicitly 
their center manifold. Due the physical differences be- 
tween solutions with standard 4D and non-standard 5D 
behaviors, we consider the cases fl\ — and fl\ 7^ in 
separated subsections. 

Summarizing, the possible late-time stable solutions 
are 

• The scalar field-matter scaling solution (P4) pro- 
vided < 7 < §,s* < -V37, /'(**) < or 
0< 7 <f,s*> V37,/'(«*)>0. 

• The scalar field-dominated solution (P5) provided 
< 7 < I, -V37 < s* < 0, /' («*) < or § < 7 < 
2,-2 < s* < 0,/'(s*) < or < 7 < |,0 < s* < 
V^7, /'(**) > Oor § < 7 <2,0< s* < 2,f («*) > 
0. 

• The scalar field-dark radiation scaling solution (P7) 
provided I < 7 < 2,s* < -2,f'(s*) < or § < 
7 < 2,s* > 2,/'(s*) > 0. 

• The de Sitter solution Pio with tt\ = provided 
/(0) is a real (finite) positive number, i.e., /(0) > 0. 

It is worth noticing that there are fractional terms with 
x 2 + y + Q, m in the denominator, in the equations (|2"Tj) - 
p6|. To finish this section we want to discuss the case 
when the denominator get close to zero and the effects 
this may have for the system. First than all, the only 
singular solution located in the hypersurface x 2 + y + 
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TABLE III: Eigenvalues for the critical points in TableU We use the notation /3± = § ( 7 - 2 ± J {2 - 7) ( -97 + 2 



Pi Eigenvalues 



Pi {-37,37,^,^,87-4,0} 

P2 {-6,6,2,0,0,6-37} 

if {-6,2,0,0,6-3 7 ,6t^6s*,=F^6s* 2 /'(s*)} 

Pa {37 - 4, fcl, _3 7 , /3+j j9 _ ) -3 7 s*/' (s*)} 

P 5 {| (s* 2 - 6) , i (s* 2 - 6) , s* 2 - 4, s* 2 - 3 7 , -s* 2 , -s* 3 /' («*)} 

P± {2, 0,0, 6 -3 7 , -6, -6s*/' («*)} 



Pt {-1,4 - 3 7 , -4, -i (l + - 15) , -I (l - ^ - 15) , -4s*/' (s*)} 

p£ {6,2,0,0,-6,6-37} 

P 9 {2, 0, -6, 6 - 3 7 , 6 - VGxs* , -V(ixs* 2 f (s*)} 

Pro {-4,0,0,-3,-3,-37} 

Pu {0,3(7-1), 3(7-1), 3 7 , 67, 67 -4} 



f2 m = is the fixed point Pu. In order to investigate the 
dynamics of (|21l) - (j2"6")) near Pu we introduce the local 
coordinates 



{x,y, fix - 1, fl a , Q m , s - s c } = 

e{x,y,n x ,n a ,n^,s} + 0(e) 2 , 

where e is a constant satisfying e<l. 



3 „ ( 3-fr 

-s c y + x 6 

2 " \l+r 



(33) 



V 



1 + r' 
3 7 r (y + O m + 2n x 



1 + r 

7r 



(2/ + ^r, 



7 -- 
where 



1 + r 
3jr ( fi m - y 

1 + r 

-V6s c 2 xf(s c ) 



-1 , 



(34) 



From the above equations it is deduced the equation r' = 
—3jr. Since 7 > 0, follows that r goes to zero (resp. to 
infinity) at an exponential rate as t -> 00 (resp. t — > 
—00). Thus, in the limit t — » —00, r/(l + r) — > 1, and 
the system (f3~4"| has the asymptotic structure 



x' = 3(7- l)ar+ y^Scy, 
?7 = 67?/, 

^' = 2(3 7 - 2)0^ + (3 7 - 4)(y + TQ, 
^' = 3(7-1)^, 

n m = 37 fn m - , 

V = -V6s c 2 xf(s c ) (35) 

The system (f55|) admits the exact solution passing by 
(x ,yo,^Ao> fi <TO!^mo,3o) at r = given by 



x(j) 



e 3(7-l)r (76 Sc y ( e 3(7+l)r _ fj + 6 J Q ( 7 + !)) 



6(7+1) 



_ „2(3 7 -2)r 



y(r) = y e' 
Sh(r] 

Kn(r) = e 3 ^ T (y (-e 37T ) +y + O^ ) 



(y + H m0 + Qao) ~ (2/0 + ^,„ )e 37T , 
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a(r) 



6s o7 ( 7 2 - 1) - s 2 J(s c ) (s c y ((7 - l)e & ~< T - 2 7 e 3 ^- 1 ) r + 7 + l) + 2^6x07(7 + 1) (e 3 ^ 1 ^ - l)) 



6 7 (7 2 -l) 



r2 CT (r) = n a0 e 



3( 7 -l)r 



(36) 



The eigenvalues of the linearization around the fixed 
points Pn, given by (|35[) . are displayed in the last row 
of table IIIII They are calculated under the hypothesis 
y <C f2 m as y, VL m — > (which is consistent with the limit 
r — > 00 as t — > — 00). From the analysis of the system 
(|35f follows that the singular point P\-y represents a Big- 
Bang singularity. It is non-hyperbolic with a 5D unstable 
manifold for 1 < 7 < 2. In this case it is a local source as 
numerical simulations suggest. For < 7 < 1 it behaves 
as a saddle point. Since the limit Q\ — > 1 corresponds to 
the domain where the brane corrections are important, 
and this regime is associated with the past dynamics, we 
conclude that the linear approximation (|34p is not valid 
as t — )■ +00. 



IV. DYNAMICS OF THE CENTER MANIFOLD 
OF THE DE SITTER SOLUTIONS 

In this section we will study the stability of the center 
manifold of the the dc Sitter solutions for fl\ = and 
for < fl\ < 1- For this purpose we can use the Center 
Manifold Theory (see the reviews [25l |26|). 



A. Case fi> 







The solution Piq with tt\ = can be a candidate to 
be a late time de Sitter attractor without 5D-corrections. 
To analyze their stability we carry out a detailed stability 
study using the Center Manifold Theory. 

Let us assume that /(0) € R \ {0}. To prepare the 
system for center manifold calculations we introduce the 
new variables 

ui = s, u 2 = A , vi = -1 + y + fl m + f2 A , v 2 = fi<r, 



v 3 = x —, v 4 = ft,-, 

V6 



(37) 



Then, we Taylor-expand the system u' 1i u' 2 , v[, v' 2l v 3 , v' 4 
in a neighborhood of the origin with error of order 0(4). 

Accordingly to the Center Manifold theorem, the local 
center manifold of the origin this vector field is given by 
the graph 

Wioc(°) = {(ui,u 2 ,v ll v 2 ,v 3 ,v i ) : Vi = fi(xi,x 2 ), 
i = 1 . . . 4, f (0) = 0, Df (0) = 0, u\ + u% < 6} (38) 

for S > an small enough real value. 

Deriving each one of the functions in (I38[) with re- 
spect r and substituting and substituting the vector field 



Ui, u' 2 , v[, v' 2 , 1*3, V4 one can obtain a system of quasi- lineal 
partial differential equations that the functions fi must 
satisfy. 

Solving this system using Taylor series up to an error 
term 0(4) we obtain 



1'3 



u 1 u 2 
~ 3 

3VE 



v 2 =0, 
Vi = 0. 



(39) 



Thus, the dynamics on the center manifold, is given by 



ui = -/(0)«? + O(4), 
u' 2 = -u\u 2 + 0{A). 



(40) 
(41) 



which is the same given by (36)- (37) in [19| with the 
identifications u\ = x±, u 2 = x 2 . 

Using the same arguments as in [l9[ we obtain that 
the origin is asymptotic stable for initial conditions in a 
vicinity of the origin whenever /(0) > 0. 

From the asymptotic stability of the origin of (|40l) - (|4"Tj) 
for the above choice of sign for /(0), follows that the 
center manifold of Pio (for il\ — 0) is locally asymptotic 
stable, and hence, the solution Pio of the system (f2~l"T) - 
(|26|) also is. 

Therefore, Piq with fix = and /(0) > corresponds 
to a late time de Sitter attractor. This result for RS- 
2 brane cosmology is in a perfect agreement with the 
standard 4-dimcnsional TGR framework. 



B. Case fix ^ 

In this section we investigate the stability of the curve 
of critical points Piq for < Cl\ < 1. This solution cor- 
responds to a de Sitter solution with 5D-corrections. 

According to the RS model this solution cannot behave 
like a late time attractor since 5D-corrections are typical 
of the high energy regime (early universe) and not to 
the low energy regime (late universe). If we can prove 
that this solution is of saddle type, this behavior can 
be correlated with a transient inflationary stage for the 
universe. In order to verify our claim we appeal to the 
Center Manifold Theory. 

Let us consider an arbitrary critical point with coordi- 
nates (x = 0, y = 1 — u c , ft\ — u c , £! CT = 0, f2 OT = 0, s = 0) 
located at P w . 

In order to prepare the system (|2i p - (|2"6"|) for the appli- 
cation of the Center Manifold Theorem we introduce the 
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coordinate change 

mi = s, u 2 = -u c (y 
vi =u c (y + Q m -4 

V3 



s(u c — l) 



V6 



- 2Q m ) - (u c - 

n A - 1) , «a = 

- X, V4 = U c U r 



- i)n A> 



(42) 



Then, we Taylor expand the system u[, u'^, v[, v' 2 , v' 3 , v'^ 
in a neighborhood of the origin with error of order 0(4). 

Accordingly to the Center Manifold theorem, the local 
center manifold of the origin for the resulting vector held 
is given by the graph: 



ir, 



(0) = 

= 1.. 



{(Ul,U2,Vl,V2,V 3 ,Vi) : V. 

.4,f(0) = 0,Df(0) 



0,ul 



fi(xi,x 2 ), 
h u\ < 5} (43) 

for (5 > an small enough real value. 

Deriving each one of the functions in (14U|) with respect 
r and substituting the vector field u' 1: u' 2 , v' 1: v 2 ,v' 3 ,v' 4 one 
can obtain a system of quasi-lineal partial differential 
equations that the functions fi must satisfy. 

Solving this system using Taylor series up to an error 
term 0(4) we obtain 



vi = -\u\(u c - l)u c (2u 2 + u c - 1), v 2 = 0, 

ui(u?(nc-l) 2 /(0)-3M 2 ) 
3V6 ' 



V A = 0. 



(44) 



Thus, the dynamics on the center manifold, is given by 

tti=«?(t» c -l)/(0) + O(4), (45) 

u' 2 = -u\ (K(3m c - 4) + l)u 2 + (u c - lfu c ) + 0(4). 

(46) 

which is the same as (41), (42) in [l9| under the variables 
re-scaling 



(1 - u c ) m 
V6 



-, u 2 



-u 2 . 



Thus, using the same arguments as in [19j . the origin of 
coordinates is locally asymptotically unstable (of saddle 
type) irrespective the sign of /(0). Henceforth, the center 
manifold of Pio is locally asymptotic unstable (saddle 
type) for /(0) ^ 0. 

The physical interpretation of this result is that there 
are not late time attractors with 5D-modifications. This 
type of corrections are characteristic of the early uni- 
verse. In this sense the solution Pio with < Q\ < 1. is 
associated to the primordial inflation. 

As in [l9| for Pio with < fl\ < 1, we assume that 
/(0) € M \ {0}. Otherwise it is required to include higher 
order terms in the Taylor expansion, increasing the nu- 
merical complexity. 



V. A TOY MODEL 

The objective of this section is to illustrate our analyt- 
ical results for the following toy model. Let us consider 
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P7 is stable 



P4 is stable 



P5 is stable 



P10 is stable 



FIG. 1: Parameter space diagram showing the sufficient con- 
ditions for the existence of late-time attractors for the input 
function f(s) = s + a. 



a simple /-deviser given by f(s) = s + a. This function 
corresponds to the potential given in the implicit form 



1 



(s) = — + In 

sa 



V(s) - V 



(47) 



or alternatively 

v{4>) = v 



(48) 



where W(z) is the special function 'ProductLog' that 
gives the principal solution for id in z = we w . V(cj>) is 
defined for (f> > 0. For a < 0, V is a monotonic decreas- 
ing function taking values in the range [0, Vb]; for a > 0, 
V is a monotonic increasing function taking values in the 
range [V ,oo). 

The function / satisfies s* — —a, f'(s) = 1, /(0) = a 
and s*f'(s*) = s* = —a. The sufficient conditions for the 
existence of late-time attractors are fulfilled easily. They 
read (see the figure [T|): 

• The scalar field-matter scaling solution (P4) is a 
late-time attractor provided 0<7<|,a< — 

• The scalar field-dominated solution (P5) is a late- 
time attractor provided < 7 < ^,—y/3j < a < 
or § < 7 < 2,-2 < a < 0. 

• The scalar field-dark radiation scaling solution (P7) 
is a late-time attractor provided | < 7 < 2, a < 
-2. 

• The de Sitter solution Pio with £l\ — is the late- 
time attractor provided a > 0. 
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FIG. 2: (a) Some orbits in the invariant set Qx = f2 CT — 0,s = 
2 of the system (|21|) -(I26 | ) for the input function f(s) — s + a 
for a = —2. (b) Projection in the plane x,y. This numerical 
elaboration shows that P4 is the local attractor and P^ 1 are 
early-time attractors for this invariant set (actually they are 
saddles in the 6D phase space.) 



FIG. 3: (a) Projection in the plane y-£l\ for the system (|21[) - 
(|26p restricted to the invariant set Q m = n CT = 0, s = 2 for the 
input function f(s) = s + a for a — —2. It is showed there the 
instability of the line y + Q\ — 1 and the stability of P4. In 
(b) are displayed some orbits projected in the plane £lcr,£l\. 
This figure illustrate that the universe evolves from a solution 
with non-standard 5D behavior to an isotropic solution (P4). 



In the figures [5] and [3] are presented some numerical in- 
tegrations for the system ([2~Tj) - (l2l)l) for the input function 
f(s) = s + a for a = —2. In this case the local attrac- 
tor is the scalar field- matter scaling solution (P4). The 
point P5 and P7 coincides in which case they are non- 
hyperbolic (with saddle- type behavior). We investigate 
the invariant set s = 2 which contains the relevant late- 
time attractor. Since the system is 6D the dynamics is 
more richer than showed in the figures. However, as we 
proved analytically, the universe evolves to an isotropic 
standard scalar field-matter scaling solution. 

In the figures 2] and [S] are presented some numerical in- 
tegrations for the system ([2"Tj) - ([2l))) for the input function 



f(s) — s + a for a = — 1. In this case the local attractor 
is the scalar field-matter scaling solution (P5). P4 and Pj 
do not exist. We investigate the invariant set s = 1 which 
contains the relevant late-time attractor. As before, the 
system is 6D the dynamics is more richer than showed in 
the figures, but, as we proved analytically, the universe 
evolves to an isotropic standard scalar field dominated 
solution. 

In the figure [6ja) are displayed some orbits in the in- 
variant set fl m = CT = 0, s = of the system (I2ip - (|2^1) 
for the input function f(s) — s + a for a = 0.1 In[5Jb) is 
presented the projection in the plane x, y. This numerical 
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FIG. 4: (a) Some orbits in the invariant set Qx = f2 CT = 0, s = 
1 of the system (|21|) -(I26 | ) for the input function f(s) — s + a 
for a = —1. (b) Projection in the plane x,y,Q m = 0. This 
numerical elaboration shows that P5 is the local attractor and 
Pg are early-time attractors for this invariant set (actually 
they are saddles in the 6D phase space). P4 and P7 do not 
exist. 



elaboration shows that P10, = is the local attractor 
and P-f are saddles. The local source in this invariant 
set is Pn. In^c) are drawn some orbits projected in the 
plane y-Q\ for fl m = suggesting that the line y+Q\ = 1 
is a local attractor in the invariant set Q m — 0. However 
they are saddles for the full dynamics. 



VI. RESULTS AND DISCUSSION 

The main results of this investigation can be summa- 
rized as follows. The singular point Pn represents a Big- 



FIG. 5: Projection in the plane y-Qx for the system (|21[) - (|26|) 
restricted to the invariant set Q m = fi CT = 0, s = 1 for the 
input function f(s) — s + a for a — — 1. It is showed there 
the instability of the line y + fl\ = 1 and the stability of P5. 



Bang singularity. It is non-hyperbolic with a 5D unsta- 
ble manifold for 1 < 7 < 2. In this case it is a local 
source as numerical simulations [5] and ^c) suggest. 
For < 7 < 1 it behaves as a saddle point . 

In the general case /(0) € K, the solution P\o with 
fix = is the late time de Sitter attractor without 
5D-corrections for /(0) > 0. To analyze their stability 
we have used the Center Manifold Theory. Using this 
technique we have obtained an analogous center mani- 
fold for the origin (up to a variable-rescaling) to the one 
in (41), (42) in [Ig|. This allow us to prove that P w 
with tt\ £ (0, 1) are not late time attractors with 5D- 
modifications since they are always saddle-like. This fact 
correlates with a transient primordial inflation. These 
results extent our results in [19| to the Bianchi models 
class. 

The critical points P 2 ± are solutions dominated by 
the kinetic energy of the scalar field, representing so- 
lutions with an standard behavior (Q,\ = 0), although 
non- hyperbolic, they behave as saddle-like points in the 
phase space. Also are the critical point P3 which are 
solutions dominated by the kinetic energy of the scalar 
field. The critical points Pg corresponding to scalar 
field-anisotropic scaling solutions (f2^ ~ f2 CT ) are non- 
hyperbolic, however they behave as saddle-like points 
since they have both nonempty stable and unstable mani- 
folds. All these cosmological solutions represent transient 
states in the evolution of the universe. 

The possible late-time stable solutions are 

• The scalar field-matter scaling solution (P4) pro- 
vided < 7 < §,s* < -V37,/'0*) < or 
0<7<t,s*> VS7,/'(s*)>0. 

• The scalar field-dominated solution (P 5 ) provided 
< 7 < f , -V37 < s* < 0, /' («*) < or § < 7 < 
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2, -2 < s* < 0, /' (s*) <0or0<7<f,0<s*< 
V^y, /' (a*) > or | < 7 < 2, < s* < 2, f (s*) > 
0. 

• The scalar field-dark radiation scaling solution (P7) 
provided | < 7 < 2, s* < -2,/'(s*) < or f < 
7<2, S *>2,/'( S *)>0. 

• The dc Sitter solution P\q with tt\ = provided 
/(0) is a real (finite) positive number, i.e., /(0) > 0. 

The main difference with respect to our analysis in [191 ] is 
that there are possible scalar field-dark radiation scaling 
late time solutions (P7) for a wide region in the parameter 
space (see figure [1] for a numerical elaboration) . 

Our results are quite general and apply to the scalar 
fields potentials presented in table I in the reference [l!| . 
In the particular case of a scalar field with potential 
V = Voe~ x ^ + A it can be proved, in a similar way as in 
(l9j . that for x < the de Sitter solution is asymptoti- 
cally stable. However, for x > the de Sitter solution is 
unstable (of saddle type). We omit the calculation here 
because it is straightforward and the result is consistent 
with our previous results in [l9l | . Instead of considering 
classical potentials discussed in the literature we have 
considered a toy model for illustrating our analytical find- 
ings. In this case we have investigated a scalar field with 



potential V(<f>) = Vq 



where W(z) is 



the special function 'ProductLog' that gives the principal 
solution for w m z — we w . This potential corresponds to 
the /-deviser f{s) — s + a. 

In this paper we are mainly interested in the late-time 
dynamics. For the above particular example the late- 
time dynamics can be summarized as follows. 

• The scalar field-matter scaling solution (P4) is a 



late-time attractor provided < 7 < |, a < —y/Ej, 

The scalar field-dominated solution (P5) is a late- 
time attractor provided < 7 < |, — \fW) < a < 
or I < 7 < 2,-2 < a < 0. 

The scalar field-dark radiation scaling solution (P7) 
a 
2. 



is a late-time attractor provided | < 7 < 2, a < 



• The de Sitter solution P10 with fl\ = is the late- 
time attractor provided a > 0. 



(c) 



VII. CONCLUSIONS 



FIG. 6: (a) Some orbits in the invariant set Q m = S1 CT = 0, s = 
of the system (I21I) - (I26I ) for the input function f(s) — s + a 
for a = 0.1. (b) Projection in the plane x,y. This numerical 
elaboration shows that P10, tt\ = is the local attractor and 
P* are saddles. The local source in this invariant set is Pn. 
(c) Projection of some orbits in the plane y-£l\ for Q m ~ 
suggesting that the line y + Q\ = 1 is a local attractor in the 
invariant set Q m = 0. However they are saddles for the full 
dynamics. 



In the present paper we have investigated the phase 
space of the Randall-Sundrum braneworlds models with 
a self-interacting scalar field trapped in a Bianchi I brane 
with arbitrary potential. 

From our numerical experiments we claim that Pn is 
associated with the Big Bang singularity type. The nu- 
merical investigations performed in this paper suggest 
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that it is in general the past attractor in the phase space 
of the Randall-Sundrum cosmological models. 

Using the center manifold theory we have obtained 
sufficient conditions for the asymptotic stability of de 
Sitter solution with standard 4D behavior. We have 
proved, using the center manifold theory, that there are 
not late time de Sitter attractors with 5D-modihcations 
since they are always saddle-like. This fact correlates 
with a transient primordial inflation. We have obtained 
sufficient conditions on the potential for the stability of 
the scalar field-matter scaling solution, the scalar field- 
dominated solution, and the scalar field-dark radiation 
scaling solution. We illustrate our analytical findings us- 
ing a simple /-deviser as a toy model. All these results 
are generalizations of our previous results obtained for 
FRW branes in 
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